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Abstract 

In this paper, we consider the solutions of the relaxed Q-tensor flow in K.^ with small 
parameter e. Eirstly, we show that the limiting map is the so called harmonic map 
flow; Secondly, we also present a new proof for the global existence of weak solution for 
the harmonic map flow in three dimensions as in |23| and [18) . where Ginzburg-Landau 
approximation approach was used. 


1 Introduction 

Liquid crystals are a state of matters that have properties between those of a conventional 
liquid and those of a solid crystal. One of the most common liquid crystal phases is the 
nematic. The nematic liquid crystals are composed of rod-like molecules with the long axes of 
neighboring molecules approximately aligned to one another. There are three different kinds 
of theories to model the nematic liquid crystals: Doi-Onsager theory, Landau-de Gennes 
theory and Ericksen-Leslie theory. The first is the molecular kinetic theory, and the later 
two are the continuum theory. In the spirit of Hilbert sixth problem, it is very important to 
explore the relationship between these theories. 

Ball-Majumdar[T] define a Landau-de Gennes type energy functional in terms of the mean- 
field Maier-Saupe energy. Majundar-ZarnescufH] consider the Oseen-Frank limit of the static 
Q-tensor model. Their results show that the predictions of the Oseen-Frank theory and the 
Landau-De Gennes theory agree away from the singularities of the limiting Oseen-Frank 
global minimizer. 

In [mis], Kuzzu-Doi and E-Zhang formally derive the Ericksen-Leslie equation from the 
Doi-Onsager equations by taking small Deborah number limit. In [201121] . Wang-Wang-Zhang 
present a rigorous derivation from Doi-Onsager theory and Landau-de Gennes theory. In [9], 
a systematical approach was proposed to derive the continuum theory from the molecular 
kinetic theory in both static and dynamic case. 
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Different with the above results on the static or local case, the goal of this work is to 
investigate the global convergence problem of the solutions of the Q-tensor flow in 
We will show that convergences weakly to the weak solution of the harmonic map flow. 


1.1 The relaxed Q-tensor flow 


In Landau-de Gennes theory, the state of the nematic liquid crystals is described by the 
macroscopic Q-tensor order parameter, which is a symmetric, traceless 3x3 matrix. Phys¬ 
ically, it can be interpreted as the second-order moment of the orientational distribution 
function /, that is 


f 1 

Q= (mm-/d)/(im. 

Js^ 3 


When Q = 0, the nematic liquid crystal is said to be isotropic. When Q has two equal 
non-zero eigenvalues, it is said to be uniaxial and Q can be written as 

Q = s(nn - Id), n G S'^. 
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When Q has three distinct eigenvalues, it is said to be biaxial and Q can be written as 

Q = s(nn — -Id) + A(n^n' — -Id), n, G 5^, n • n' = 0. 
o o 

The general Landau-de Gennes energy functional takes the form 

J^cg{Q,^Q) = [ |-VrQ2-brQ3 + ^(trQ2)2 

Jr3 _ 3 ^ 4 _^ 

/B:bulk energy 

■77 (-^ll^QI T L 2 QijjQik,k T L^Qij kQik,j T I-“iQijQkl,iQkl,j^ "ifix, 

Z s-J 

/B:elastic energy 




here a, b, c are material-dependent and temperature-dependent nonnegative constant and 
Li{i = 1, 2, 3,4) are material dependent elastic constants. We refer to [H |T3] for more details. 

There are several dynamic Q-tensor models to describe the flow of the nematic liquid 
crystal, which are either derived from the molecular kinetic theory for the rigid rods by 
various closure approximation such as pmE], or directly derived by variational method 
such as Beris-Edwards model [3] and Qian-Sheng’s model HZ!. 

In [ 21 ], the authors consider the following Beris-Edwards model 


(Q) 




+ V • Vv = -Vp -k V • (o-^ -h cr“ -k cr'^), 


dt 

V • V = 0, 

-kv-VQ-kQ'D — D-(5 = —H-\- Sq{D). 


( 1 . 1 ) 


Here T is a collective rotational diffusion constant, and 


D = l(Vv + (Vv)^), rj = i(Vv-(Vv)^). 
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Moreover, , cj“ and are symmetric viscous stress, antisymmetric viscous stress and 
distortion stress respectively defined by 


a* = 77 D - 5q(H), a“ = Q.H-H-Q, 
where 77 > 0 is the viscous coefficient, H is the molecular field given by 


H(Q) 




( 1 . 2 ) 


and S'q(M) is defined by 

5q(M) = e (^M • (Q + il) + (Q + il) • M - 2(Q + il)Q : m) 


for symmetric and traceless matrix M, where ^ is a constant depending on the molecular 
details of a given liquid crystals. The well-posedness results of the Q-tensor model are studied 

in [T5l [16] . 

Wang-Zhang-Zhang [2T] justify the limit from Beris-Edwards system with a small param¬ 
eter e to the Ericksen-Leslie system before the first singularity time of the limit system. The 
limit behavior of the solution after the singularity remains unknown. In this paper, we are 
interested in the global convergence from the Q-tensor flow to the harmonic map flow in 
Let us begin with the simplest form L 2 = L 3 = L 4 = 0, i.e., 

^ = -|;W) + ^AQ, (1.3) 

where 

J{Q) ■■= = -aQ - bQ‘^ + clQpQ + IftlQpI. 


Let 6 e 5^ is a constant vector, no : —)• 5^ such that no — 6 £ LI®+^(M^)(s > 0). 

Moreover, Qo{x) = s+(no(x) (g) no(a:) — ^). We consider the following relaxed Q-tensor 
equations with a small parameter e: 


(Qe) 


dQe _ (® “ c\Qe\‘^)Qt + hQl - LiAQ^ 

< 'W ~ ef ^ f 

^ Q6(-,0) = Qo{x), 


(1.4) 


which has a unique strong solution Qeit, x) satisfying Qe{t, x)—Qo{x) £ L^^([0, 00 ), 
We will study the global convergence of as e tends to zero. 


1.2 Main result 

The initial data Qq of the Q-tensor flow equations ()1.4I) lies in a special space, which contains 
the minimizers of the bulk energy fsiQ)- To begin this, we introduce some notations and 
known results. 
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Let denote the set of real 3x3 symmetric matrices and Qq C denote the 

space of Q-tensors defined by 

Qo :={Q&M^y^,Qu = 0}, 

where we have used the Einstein summation convention. The matrix norm is defined as 

IQI := \/trQ2 = 

We also write 


|vq|2 = ^^Q,j^aQ^y 


The bulk energy density can be written as 

fB{Q) = -\\Q?-\tT{Q^) + l\Q\^. (1.5) 

One can verify that Jb is bounded from below (for example, see |14l Proposition 8]), thus Jb 
has the corresponding non-negative bulk energy density Jb defined by 

fsiQ) = fsiQ) - min /i?((5). (1.6) 

Q^Qo 

In [141 Proposition 8], it was proved that Jb attains its minimum on the uniaxial Q—tensors 
with constant order parameter as shown below 


/b(Q) = 0 Q G M where 


Af = \ Q G Qo, Q = s+ ( n (g) n - ) , n G 5^ 


with 


■S+ = 


b + Vb'^ + 2Aac 
4c 


(1.7) 


( 1 . 8 ) 


For a matrix Q G Af, we use TgAf to denote the tangent space to AA at Q m Qo) and (TqAA)^^ 
to denote the orthogonal complement of TgAf to Qo- 

Let {A, B} = AB -|- BA for A,Bg It was described in [221 Lemma 2] that 


TqM = \^QGMll^-.^s+Q = {Q,Q}Y 

= {n (g) h + h (g) n : h, h G TnS^}, (1.9) 

and 

{TqMA = {q^gQo:Q^Q = QQ^}. ( 1 . 10 ) 

We will specifically describe the orthogonal basis for the tangent and normal space in the 
following Lemma 1221 It is obvious from (jl.lOp that J{Q) G {TgAf)^^ for Q G Af, and we 
will show in Lemma 12.41 that for the approximating near Af, we still have 

J{Q,) G 
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where denotes the projection operator on M. 

Let z = (x,t) denote points in x M+. For zq = {xo,to), R > 0, let 

Pr{zo) = {z = {x,t)\\x - xo\ < R,\t- to| < R^}, 
Sr{zo) = {z = {x,t)\t = to- R^j, 
Tr(zo) = {z = ix,t)\to - AR^ <t <to- R^}. 


Denote the scaled fundamental solution to the heat equation 


where 


Then we have 


Gzaiz) = Gf T ^ r , N f ) , 


G.o(z) = 


(47r(to - t))y^ 


e , t < to- 




r X - xq ^ 
2Li t — to 


Also we write Pr{0) = Pr, Tr{0) = Tr, and Go{z) = G{z). 

Similar to the harmonic map flow, the limiting Q-tensor flow takes as follows 


dtQ - ^AQ + Xix,t)'yj^{Q) = 0, Q(x,t)|i=o = Qo(x), (1-11) 

where A(x,t) is a function of and 'y^f{Q) is unit normal vector to (TqAA)q^ at Q. 

Definition 1.1 A Q—tensor Q{x,t) : x M_|_ M is called a weak solution to if 

Q{x,t)\t=o = Qo{x) a.e., dtQ,VQ G x M+), and there holds 



{dtQ : 4> + : V4> + XjjviQ) ■ 4>)dxdt = 0, 


for all (j) G x M+,M3 x3)_ 


Our main theorem is stated as follows. 


Theorem 1.2 Let Qe satisfy the equations of the relaxed Q-tensor flow equations with 
the data Qo{x) = s+ (no <8) no — ^Id) ^ N as in Then 

1. there exists a subsequence of e (also denoted e) such that 


Qe ^ Q 

and Q solves weakly the equations 


n I 


n- Id 

3 


eM, 
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2. the director field n weakly solves 

Sin - An = |Vn|^n, n(x, i)|i=o = no(x). (1-12) 

3. n is regular on a dense open set Oq C x M_|_, whose complement S has locally finite 
3-dimensional Hausdorff-measure (with respect to the parabolic metric). 

Remark 1.3 Compared with I21A 121^ . the above theorem makes it reasonable that the global 
weak convergence from the Q-tensor flow to the simplified Oseen-Frank map flow (i.e, har¬ 
monic map flow). Different from Ginzburg-Landau approximation used in and f7^ . we 
consider the Q-tensor approximation, and the difficulty is that the properties of the limit 
manifold M is unclear as stated in f2E/ . In the next section, we give a careful study for the 
geometry of M ( see Lemma 


2 Technical lemmas and interior regularity estimates 


In this section, we will introduce the properties of Q-tensor matrix, the tangent space, the 
normal space and the equivalence of the bulk energy. Using these estimates and exploring 
monotonicity inequalities as in [23], we can obtain the interior regularity criteria of Q-tensor 
equations. 

First of all, for the matrix of 3 x 3, we have the following properties. 


Lemma 2.1 Let A,B he matrices o/3 x 3. 
(i) If A is symmetric, then 


A: B = A: B, 
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where B is the symmetrization for B. 

(a)If A is antisymmetric, then 

D _ dT 

A: B = A: B, B = - 

’ 2 

where B is the antisymmetrization for B. 

(Hi) If A is symmetric and B is antisymmetric, then 

A:B = Q. 

For Q G Af, it is easy to verify that the orthogonal basis of TqAI and is as 

follows. 

Lemma 2.2 Let Q = s+ (ns (g) ns — \ld) G M, and ni,n 2 he unit perpendicular vectors in 
Vhg = {n-*- G : n-*- • ns = 0}. Then it holds that 


1 . 


TqM 


Span 


I ^(ns (g) n2 -I- n2 (g) ns), ^(ns (g) ni -|- ni (g) ns) 






( 2 . 1 ) 
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2 . 




Span 

Vq 


1 


(n2 (g) ni + ni 


N 1 . 

n 2 ), ^(ni (g) ni - n 2 (g) n 2 ), 
V2 


-ni 


1 

ni + -n.2 


n2 


Id\ 

"tJ 


( 2 . 2 ) 


5. Moreover, Qq = TqM ® {TqM)^^. 

Proof: The tangent space at s+(n 3 (g) na — i/d) of Af ( 12 . 11 ) is a direct result from [211 (2.3)]. 
The others can be deduced by direct computations. □ 


Lemma 2.3 For Q G Qo; there exists eo > 0 such that if dist{Q,A/) < cq, then 

JiQ) e (2-3) 


Proof: Denote the eigenvectors of Q{x,t) by 111(3:, t),n2(x,t),n3(x,t) corresponding to its 
eigenvalues Ai(3:,t), X2{x,t), Xs(x,t) = —Xi{x,t) — X2{x,t). Then we have 

Q{x, t) = Aini (g) ni + A2n2 (g) n2 + Aana (g) 113, ( 2 . 4 ) 

especially, 

/ = ni (g) ni + n2 (g) n2 + 113 (g) 113. ( 2 . 5 ) 

Choose eo small enough such that dist{Q,M) < eo, then 

dist{Q,Af)‘^ = ^Ai H—+ ^A2 H—+ ^Ai + A2 + 2 “^^ > ( 2 - 6 ) 

furthermore, 

dist{Q,J\ff = \Q- TTAf{Q)\‘^, ’n'AfiQ) = 'S+(n3 (g) no - y) 

where T^'jviQ) is unique and depends continuously on Q. See | 22 l Lemma 8] for more details. 
It is easy to see the projection of Q on is 0 . Using () 2 . 5 p . we have 


Q 


Aini (g) ni + A2ni (g) 112 - (Ai + A2)n3 (g) 113 
( 2 Ai + A2) f ni (g) ni - y j + (Ai + 2A2) [^2 ® ^2 



113 


= Aini (g) ni + A|n 2 (g) n 2 + (Ai + A2)^n3 
= — A 2 ( 2 Ai + A 2 ) ^ni (g) ni —— Ai(Ai + 2 A 2 ) 

(Af + Af + A 3 ) 


n2 


Id 

n2- 

3 


+- 


-Id, 
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which together with (I2.4|) yields that 


IQI^ 


Id 


Q - — —A2(2Ai + A 2 ) { ni (g) ni —^ ) — Ai(Ai + 2 A 2 ) ( n2 (g) n2 —— j . 


Id\ 


Hence, we get 


J{Q) = (a-c|g|2)g + 6 




Id 


Id 


Id^ 


= (a - c\Qf - 6 A 2 )( 2 Ai + A2)(ni (g) ni - y) + (a - c\Q\^ - 6Ai)(Ai + 2A2)(n2 (g) n2 - y) 
= (a - c|gp - 6A3) (Ai - A2)e2 

H—^ ((o — c|gp — 6A2)(2 Ai + A2) + (a — c|gp — 6Ai)(Ai + 2A2)) 63, 

V6 


where we have used the orthogonal basis 62 and 63 in the normal space {Titj^{Q){-^))qq 
(| 2 . 2 I) ). and 


see 


1 / N 1 

62 = ^(ni (g) ni - n2 (g) n2), 63 = V6(-ni (g) ni + -n2 (g) n2 - y). 

Thus, (j2.3p is an immediate result. □ 

In fact, the nonnegative bulk energy fsiQ) is equivalent to the distance from Q to M, 
which is stated as the following Lemma. 

Lemma 2.4 There exists cq > 0 such that if dist{Q,J\f) < cq, then 

^dist{Q,Aff < /s(g) < C{dist{QM))\ (2.7) 

where C is independent of Q, but depends on a, b, c. 

Proof: Assume that the eigenvalues of Q are x,y, —x — y. If dist{Q,M) < cq and cq is small 
enough, similar to (12.61) . we have 

dist{Q,J\f)‘^ = (x + y+ (y + + [x + y + 2y< ei- (2-8) 

Let 

{x + y + {y+ y + {x + y+ = G{x, y). (2.9) 

On the other hand, for the nonnegative bulk energy, we have 

fsiQ) = -^|gp-^trg3 +jigi"^- min /s(g) 

2 3 4 QgSo 

= -a{x^ + + xy) + b{x^y + xy'^) + c(x^ + y^ + xy)'^ - min fsiQ) 

QgQo 

= H{x,y), (2.10) 
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where H : —)• M is the 2-dimensional function as G{x,y). Note that H{x,y) = 0 only at 
three pairs {x,y) namely (2^,—^) and (—^,2^), c.f.[TU Lemma 5]. This 

fact together with (12.9p and ()2.10p gives 


( 2 . 11 ) 


' ■S+ 


dH / 

s+ 

. 3 ’ 

3 ) 

dx V 

3 ’ 

f_£± 

S+'' 

) = —1 

r_£± 

V 3 ^ 

’ 3 / 

' dx 

1 3 



' s+ 


) dy \ 

. 3 ’ 

3 ) 


f_£± 

_£±^ 

) dy 

1 3 ’ 

3 / 


Careful computations show that 

dH 
dx 
d‘^H 


dxdy 

52 if 

5x2 


= (2x + y)[-a + hy + 2c{x‘^ + + xy )], 

= [-a + hy + 2 c(x2 + y"^ + xy)] + (2x + y)[b + 2c{2y + x)], 


= 2[-a + by + 2c{x‘^ + 2/^ + xy)] + 2c(2x -h y)^- 


Noting 2cs^ — 6s+ — 3a = 0, we have 


52 iL 


5x2 


i 


d'^H , 



V 3 ’ 

3 / 

dy"^ ' 

. 3 ’ 

3 / 


d‘^H 

dxdy 


(—= —bs-\. + 2cs\ = 3a, 


3 ’ 3 / 


52G , 52G 52^ 

= 4 = -^, = 


5x2 


5y2 ’ dxdy 


Then, for (x,y) 7 ^ by the above computations and (12.lip , we have 

H{x,y) Hi{x,y) + RH{x,y) 


G{x,y) 


Gi{x,y) 


where 




3 / 

le 
we get 


Gi{x,y) = i{x + :^] +4(i + ^)(9 + :^)+4(» + ;7 


•5+ 


'S+\ 




3 / 

+ 

3 


and Rh is the remainder in the Taylor expansions of H{x,y) at — 

sufficiently small in 

bs 

'2 


( 2 . 12 ) 


Hiix, y) = {bs+ + 3a) (^(x y)^ + iv + y)^) + (^x + y) (v + y) > (2-13) 


{Rnix, y)| < ^ ((x -h y)^ + (y + y)^) • 


(2.14) 
Thus for €1 

(2.15) 


Summing up the inequalities (j2.12p - (j2.15p . we conclude that 

bs+ H{x,y) 3bs+ 


The proof is completed. 


□ 
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Now we consider the evolution of the energy, and we will follow the same line as in |23] . 
First of all, we define the energy density by 

e.(Q, VQ) = ^fsiQ) + ^|VQP. (2.16) 

For the equations of (II.4p with initial data Qq G J\f, there exists a global weak solution 
Qe{x,t) (denoted by Q for simplicity), see [16]. The solution Q is regular indeed by usual 
energy estimates, and we have the following basic estimates. 


Lemma 2.5 (Energy Inequality) Suppose that Q{x,t) solves {1-4^ with initial data Qq G M, 
then it holds that 


sup 

t>o 


[ [ \dtQ\‘^dxdt + [ e^{Q{-,t),VQ{-,t))dsi 
Jo Jr3 Jr3 


< ^ / iVQoPdx. 

2 i Im 3 


(2.17) 


Proof: Multiplying (11.41) by Qt and integration by parts yield that 

[ \dtQ\^ = -, 

JiR3 el 


: dtQdx + ^ / AQ : dtQdx. 


eF J-g3 SQ F 

Noting that : dtQ = OtfsiQ) and /b(Qo) = 0, the lemma follows. □ 

The following parabolic maximal principle lemma is similar to Proposition 3 in where 
the elliptic case was considered. We omitted the proof. 

Lemma 2.6 (Maximal Principle)Suppose that Q{x,t) solves ^-4^ with initial data Qq G 
then it holds that 


\Q\ < Y 2®+. 

For the case of L 2 = L 3 = L 4 = 0, we also have the monotonicity properties of the level 
energy as in [23] . 

Lemma 2.7 Suppose that Q{x,t) solves ^1.4^ with initial data Qq G M. For any point 
zq = {xo,to) G X M+, the functions 


HR,Q,e) = ^R^ [ + M91]G,,dx, 

1 Jsh{zo) 2 e 


^(i2,Q,e) = i [ [h\^Q\^ + M91]G,^dxdt 

^ JTn(za) 2 e 


(2.18) 

(2.19) 


Mji{zo) 

are non-decreasing for 0 < 72 < 

Proof: We first note that Q{x + xo,t + to) satishes the equations of (jl.4p with initial data 
Q{-,—to) = Qo{x). Thus, we may assume that zq = (0,0). By scale invariance QR{x,t) = 
Q{Rx, R^t) satisfying (11.41) with constant eji = we have 


$(i2,Q,e) = ^R^ [ + i^]Gdx 

^ 2 e 




2 , /b(< 3 -r) 


+ 


e/722 


]Gdx^^{l,QR,e/R^). 
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It suffices to consider the case of i? = 1. Direct computations and the equations (11.41) show 
that 


J Ili{-AQ) : {x-VQ + 2tdtQ) + ^^ : {x ■ VQ + 2tdtQ)j G{x,-l)dx 


1 

1 


r 


'Si 


Li{x • VQ + 2tdtQ) '■ S/kQ'^kGdx + 


er 


'Si 


/bG(x, -l)dx 


1 


f -dtQ : {x ■ VQ + 2tdtQ)Gdx + 

Jsi f 

[ fBG{x, -l)dx 
JSi 

{2tdtQ + X ■ VQ)'^ Gdx +-^ f 

Js 


>Si 


DG 


(x • VQ) : (x • VQ + 2tdtQ)dx 


+ 


er 


/Si 


2\t 


er 


fsGdx > 0, 


Si 


which implies the first inequality ()2.18p for R < 

For Q < R -^/^Q/2, W6 consider the term (^, e) with t Jt — R\l R") then 


^[R,Q,e) = i; [ [^iVQp + M^lGdxdt 

JTr ^ e 

[— |VQ|2 + M91jGdxdt 


1 

f 

1 


>-R 


r y-4ii;2 _ „ 

j-2R 

= 2 r~^^{r,Q,e)dr 

Jr 

|^2R 

= 2 


Jr ^[r,Q,€) 

where we used the monotonicity inequality (12.181) . 


□ 


Remark 2.8 The above lemma indicates that the monotonic radius of $ and 'I' depends on 
to, which is reasonable since we have no definition for t < 0. Similarly, if we consider the 
Q-tensor flow in x {—ARq,Rq), then ^{R,Q,e) is nondecreasing for 0 < ii < 2Ro o,nd 
'^{R, Q, e) is nondecreasing for 0 < R < Rq. 

We have the following Bochner-type inequality. 

Lemma 2.9 Suppose that Q{x,t) solves with initial data Qo € M. There exist eo > 0 
and a constant G > 0, independent of e, such that 


{dt - A)eflQ,VQ){x,t) < GeflQ,VQ)\x,t), 


( 2 . 20 ) 


provided that there exists a ball Bpj^„.-^{x) with p{x) > 0 such t/iat supj^g^^^^^^^,) dist{Q{y,t),M) < 

eo- 


11 











Proof: Direct calculation shows that 



2r 


|VQp + 


/b(Q) \ 

) 


Li 


fB 


Li 


Li 


= 9, ( ^|VQ|2 + ^ ) - ^A(|VQ|2) - ^AifsiQ)) 


2r 


er 


€r2 


: dtVQ + ^^:dtQ-^ (2VQ : AVQ + 2\V^Q\^) 


Li 6fB 


L, 


SfB 


er2(ig'^^ er2'^^5Q 


e2r2 


SfB 


6Q 


2Li 6fB. 
eV^ ^SQ 


2r2 

: VQ 

: VQ-3|V2q|2. 


( 2 . 21 ) 


It suffices to estimate the second term of the above equality. 

Denote the eigenvectors of Q{x, t) by ni(x, t), n 2 (a;, t), n3(3:, t) corresponding to Ai(3:, t), X 2 {x, t), 
X 3 {x,t) = —Xi{x,t) — X 2 {x,t). Then 

T^AfiQ) = 'S+(n3(a:,t) <8) n3(x,t) - ^/d), 

which is a minimizer of the bulk energy Jb- By the Taylor expansion of — near 'Kj^iO), 

Oy^ijOL^rnn ^ ' 

we get 


ir^l ! ir^w 

[Q[x,t)) = — [T^AfiQ)) 


dQijdQ 

mn dQijdQ 

mn 

d^h 


+ 


dQijdQ 

mn dQ 


i'^Af{Q))iQpq{^ 7Tj^[Q)pq) 0 {\Qpq{x,t) 'Xj\f{Q)pq\ ), 


pq 


where we have used the formula (11.61) of fB{Q) ■ 

Using the convex property of /s(g) at Trj\f{Q) and the maximum of Q in Lemma l2.6[ we 
have 

6 f R 9 ^ f R 

-V(-^):Vg = - QQ, ,QQ — {Q{x,t))Qmn,k{^,t)Qij,k{^,t) 

< C\Q{x,t)-T:M{Q)?+ C\VQ\^. 

Moreover, due to Lemma 12.41 we get 


-V(^):VQ < C'/B(Q(x,t)) + C|Vg|^ 


Combining (j2.2ip and (j2.22|) . we obtain 

{dt - —A.)ee{Q) + 
The proof of the lemma is completed. 


SfB 


SQ 


< Ce,(Qf. 


( 2 . 22 ) 

(2.23) 

□ 
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We consider the local uniform regularity property of the solution which follows from 
the monotonicity Lemma l2.7l and Schoen’s trick, c.f. |23l Theorem 5.1] or |19l Theorem 2.2], 


Lemma 2.10 Suppose that Q^{x,t) solves o in X [—4iig,iiQ], and there exist positive 
constants Cq and ei, such that when e < eg, for some 0 < R < Rq the following inequality 
holds 


^{R) = ^{RMe,e) = [ ( 


'Tr 


2T' 


eT 


Gdxdt < ei, 


then 


(2.24) 


I ^l|V7^ |2 , fsiQe 
sup ^|VQ,| +- 

Psr\ ^ 


< Ci5R) 


-2 


(2.25) 


where the constant (5 > 0 depends only on ee{Qo,VQo) and min{i?, 1}. 

Proof: We follow the same line as in |23l Theorem 5.1]. Let ri = 5R, 6 G (0,1/2) to be 
determined later. For r,a ^ (0, ri), r + a < ri, and any zq = (xq, Iq) £ Pr, we have 

I{a,zo) = cr~^ e^{Q^,VQe)dxdt <C e^{Qe,VQe)G(^^^^t^^ 2 P)dxdt 

JPa{zo) JPa{zo) 

— C ^e{Qej^Qe)G(^xo,to+2p) dxdt. 

JTa{to+ 2 p) 

Moreover, apply Remark 12.81 choose S small enough, and take to + 2cr^ — 4R^ = —R? and 
to + 2 it^ — 4i?2 = —4R?, then we deduce that 


I[a,zo) < Cmin 


tQ+2p—Rf i'to+2p—R: 


PiQ-rZc Jxj r 

, ^tn+2(T^—4i2? Jtr 


to+2(T^ —J to +2(j^—4i?: 


{e,{Qe,VQe)G^xo,to+2P)) dxdt 


< 


G / S-f:{Qf:,^Qf)G(^xO,to+2P) dxdt. 
Jtr 


Direct calculation shows that for given €2 > 0, if h > 0 is small enough, then we have 
G(xQ,to+ 2 P){xp) < Cexp G(x,t) 


< 


CG(x,t), if \x\ < 

0 

if |xl > ^ 
0 


GR-^ exp(-C<5- 

< CG(x,t)+ Ci?-2exp(-lnR-C<5^) 

< GG{x,t) + e 2 R~‘^, 

which holds on Tr, and here G is independent of 5 and R. Select 5 ~ (| lni?| + | lne 2 |)“^/^ 
for a small R and independent of ii if ii > 1. Thus, it follows that 


a ^ I ee{Qe,'VQe)dxdt < G4t{R) + Ge2ee{Qo,'VQo) < G{€i + e2ee{Qo,'VQo)). (2.26) 


! Per [zo) 
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For simplicity, we let e^{Q^,'VQ^) = e^{Q^). Since is regular, there exists G (Oj^) 
such that 


(ri - fJe)^ supee((5e) = max (ri - cr)^ supee(Qe). (2.27) 

Pa, 0<<7<ri 

Also, there exists a point (x£,A) G Pa, such that 

supe,(Q,) = e,(Q,)(a;e, A) = e,. (2.28) 

Pa, 

Set pe = ^(ri — cTg). Then it follows from (j2.27p and (j2.28p that 

sup CeiQe) < sup C^iQe) < 

^Pe(^€)ie) ^o'e+Pe 


Denote 


Qe{x,t) 



p x^, -h te), 


which solves the equation (I1.4D in P^e with e = ee^ and = y/elPt- Moreover, satisfies 

e,-(Qe)(0,0) = 1, supee((5J<4. 

Pr, 

If re > 1, dist{Qe,M) convergence uniformly to 0 on Pi, and there exists Cq, such that 
dist{Qe,M) < eo for e < Cq. Thus, Lemma [2T] implies that 

{dt - A)ei{Qe) < ciCeiOe), on Pi. 


Moser’s Harnack inequality shows that 

1 = ee{Qe)i0,0) < C [ ei{Qe)dxdt, 

JPi 

while, ()2.26p tells us 

/ ei{Qe)dxdt = {^/e'ef / e{Qe)dxdt < c(ei + e2eo(Qo)), 

J Pi Jp 1 (x,,t,) 

\/e7 

which leads to a contradiction if ei and €2 are suitably small. 

Hence, we may assume that < 1. Then 

1 = ee((5e)(0,0) < Cr~^ / ee{Qe)dxdt = Cr~^p~^ / ee{Qe)dxdt, 

J Pr, Pp, 

and using (|2.26l) . we get 

plee = re <C, 


then 


max {re — u)^ supe((5e) < < C, 

0<(T<r, p 

X <T 

which implies the required result by choosing a = = bP. 


□ 
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3 Proof of the main theorem and the equation of n 

By Lemma [2. 5 1 we know that for given smooth data Qq ^ M with VQo £ there 

exist a subsequence of (also denoted by Q^) and a function Q{x,t), such that as e —)■ 0, 
we have 


VQ, weakly* in L°°([0, oo); 

dtQe dtQ weakly in L^(M^ x M+), 

Qe ^ Q weakly in x 

MQe) ^ 0, in X M+), (3.29) 

which yield that 

dtQ e l 2 (m 3 X M+),VQ G L'^([0 ,oo),l 2 (m 3))^ ( 3 . 3 O) 

and hence also —)> Q ^..e. on x M+. Also there is a lifting map n G iL^(]R^) such that 


Q = s+ 



G Af. 


(3.31) 


3.1 Proof of Theorem II.2t the limit Q-tensor equations 

We follow the standard arguments as in [23] or [Ij. Define the singular set by 

^ = ni?>o{^;o G X M+l liminf / ee{Qe)dxdt > ei}. 

^^0 Jtr(zo) 

Then as in |23l Theorem 6.1](see also [I]), one can show that S is closed and has locally finite 
3-dimensional Hausdorff-measure with respect to the parabolic metric. 

For zq ^ S, there exists a > 0, and a subsequence of e, which is still denoted by e, 
such that 


e^{Q^)Gz^dxdt < ei. 

(^o) 

It follows from Lemma 12.101 that 

|VQ,1,M^ < C 
e 

hold uniformly in a uniform neighborhood D of zq. Let Q be the weak limit of (13.2911 . Then 
there exists a subsequence which we denote as Qe again, such that 


Qe^Q, in 

VQ, ^ VQ, weakly* in L^^(D). 


Note that fsiQe) < de is a polynomial of Qe, then the convergence shows fsiQ) = 0, i.e., 
Q G M. Also dist^{Qe,Af) < IQe — Qp < cq for e small enough. Then ()2.23p shows that 



GiQe) + 


1 


sJb 

6Q 


2 

{Qe) < C 


on D, 


(3.32) 
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which implies that 




6fB 


SQ 


{Q^)(f)dxdt < C{(j)), 


for all (f) G C“(n). Moreover, -^-^{Qe) = -^J{Qe) is uniformly bounded in also 

{dt — ^/^)Qe is uniformly bounded in and similar arguments hold for dtQe and 

Then we may assume that 


{dt - i^A)Qe {dt - ^A)Q, weakly in 




5h 

6Q 


m 


A, weakly in 


(3.33) 


The convergence (13.321) shows that Q{x,t) G C{Q,J\f), then the lifting n in (|3.31l) satisfies 
n G C{Q,S‘^). Lemma [231 tells us that 

s+{nl(g>nl - y) = 

where rig is the main eigenvector of Q^. 

Now we want to prove that the limit Q satisfies the equations of the Q-tensor flow (II.lip . 
By the definition of the matrix norm, the vector rig can be estimated as follows: 

lIVngllLtx. < c||V7r(Qe)||Loo < c||V(5e||oo, 

As in [21 Theorem 2], we can assume rig —)■ n uniformly on fl' C 11. 

let n| and n| be unit perpendicular vectors in which also continuously depend on 

Qe- Then the following three vectors are the basis of ■ 


ei ix,t) = (n| (g) + n"g 


n. 


2 ) ) 


el{x,t) = 


V2 


n 


2J 5 


,, . /X /1 f f 1 f f Id\ Vs, ^ Id 

egPx, t) = Vs -rig ( 8 ) rig + -n2 (g) n2 - — = —— (rig (g) rig - — 


Vs, 


3 J 


Then 


J{Qt) = f{{x,t)e\{x,t) + fl^{x,t)e\{x,t) + fl{x,t)e\{x,t). 
Note that (|3.33p also shows that 


eT 


X\{x,t), weakly in LLm, 


for i = 1,2,3. Due to Qe —> Q Cfoc(^) uniformly, we can assume nl —5- 112 and nf —ni, 
where 112 and ni are perpendicular on TnS^. Thus, 

el{x, t) ei{x, t), in i = 1, 2, 3, 


IS 

















and 


In' 


J{Qe) 

eT 


4>dxdt = [ {el — ei) ■ 4>dxdt + f — e 2 ) ■ (f>dxdt 

Jw er er 


+ 


L 


er 

eV 


(63 — 63) • (pdxdt —)• 0 


for any vector filed (j) € such that (j){z) € TqM a.e. on Q, and any il' d n, which and 

(13.3,111 yield that 


{dt - ^A)Q ± TqM, 


a.e. in fl. 


Then there exists a unit normal vector hied 7 a/'(Q) £ (^Q-^)qo Q a scalar function 

A G L‘f^^{Q) such that 


dtQ - ^AQ + XjM) = 0 (3.34) 

a.e. on n and in the sense of distribution. 

Standard arguments show that Q also weakly solves (|3.34l) on x (for example, see 
123 ]), and we omitted the details. 


3.2 Proof of Theorem II.2t the harmonic map flow 

Note that n G x M, S^) is a lifting such that Q = s+(n(8) n — -y) and |(5| is a constant. 

For 4> G X M+,M^), 


0 = 


{{(j) ■ n)dtQ : Q TVQ ■ n)Q)) dxdt + f f {4> ■ n)X'jf/{Q) : Qdxdt 

4 4m3 

\VQ\'^{(p ■ n)dxdt + f j {(j) ■'n)X'j_\f{Q) : Qdxdt. (3.35) 

Jm.^ Jm.+ 


We use (f)"' to denote (n • ()))n. Then 


XY ® {cj) - (jT) + {(j) - cjT) ® TO. ^ TqM. 


(3.36) 
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On the other hand, by (13.341) . (I3.36|) and Lemma [XT] we derive 


0 = 


{dtQ - AQ + X'yAfiQ)) : (n (g) </> + 0 (g) n) 


Q : dt{n 0 (j) + 4> ^ n)dxdt + 


VQ : V(n(g)(^ + (/)(g) n)dxdt 


+ / / X^f^f{Q) : {n0 4> + n)dxdt 

Jr^Jr+ 

= —2 / Q : {dtn (g) (/) + n (g) dtcp) dxdt + 2 / VQ : (Vn (g) 0 + n (g) V4>) dxdt 

Jr3Jr+ iR3 Jr+ 

+ f f Xjj\f{Q) : (n (g) (/)”' + (/)”' (g) n)dxdt 
iR3 Jr+ 

= —2s+ I I -n ■ dt(j) —-dtn ■ (pdxdt + 2s+ f f \Vn\‘^n ■ (p+ Vn ■ V(j)dxdt 

Jr3 Jr+ 3 3 Jr3 

+2 / / (p ■ nXjj\f{Q) : n (g) ndxdt. 

Jr3 Jr+ 

= —2s+ I I -n ■ dt(p —-dtn ■ (pdxdt + 2s+ f f |Vn|^n • (/> + Vn • V(/>(ixdt 

Jr3 Jr+ 3 3 Jr3 

+2 f [ <p ■ n—'yj\f{Q) : Qdxdt. 

JR^ Jr+ 'S+ 

= 2s+ / / dtn ■ (pdxdt — 2s+ / / |Vn|^n • ijidxdt + 2s+ / / Vn-V(pdxdt, 

Jr3 Jr , Jr3 Jr , Jr3 Jr , 


where in the last step we used (13.351) . 

Thus the proof of the Theorem 11.21 is completed. 


□ 
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